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We study the quantum dynamics of conversion of composite bosons into fermionic fragment species 
with increasing densities of bound fermion pairs using the open quantum system approach. The 
Hilbert space of A^-state- functions is decomposed into a composite boson subspace and an orthogonal 
fragment subspace of quasi-free fermions that enlarges as the composite boson constituents devi- 
ate from ideal boson commutation relations. The tunneling dynamics of coupled composite boson 
states in confined systems is examined, and the appearance of exceptional points under experimen- 
tally testable conditions (densities, lattice temperatures) is highlighted. The theory is extended to 

fT^ ■ examine the energy transfer between macroscopically coherent systems such as multichromophoric 

T-H ' macromolecules (MCMMs) in photosynthetic light harvesting complexes. 
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<^ ■ I. INTRODUCTION 

0\ ' 

Composite quasi-particle systems such as excitons (coherent superpositions of electrons and holes) display phase- 

i""!, space filling eff'ects when the mean separation between particles becomes comparable to the exciton Bohr radius. This 

"^ ' arises due to the Pauli blocking [l|, y] of scattered fermionic particles which constitute "cobosons" [3, |j] or composite 

II I bosons. The difference between composite bosons and elementary bosons can be seen for instance, in the invariance 

"^ ■ of Pauh scattering during exchanges between specific fermionic species that can be correlated via several routes Q ■ 

jrt ' Such interactions are non-existent in a collection of ideal bosons with negligible overlap, spin and internal structure. 

^ , Unlike the conventional commutation relations satisfied by ideal boson operators of bound fermion pairs, composite 

O^' bosons systems obey a series of commutation relations [a, Q that are refiective of the underlying fermionic structure of 

'""', the constituent particles. The cxcitonic polariton 0,13 is one such example of a composite boson system in which the 

I • influence of the composite-particle effect on many-body physics can be studied. When the dimensions of quantum dot 

^ , excitonic polaritons are decreased, anti-bunching features appear. The system undergoes a crossover from bosonic to 

^sD ' fermionic display in optical features, characterized by a shift from the Rabi doublet to a MoUow triplet in the optical 

OO \ spectrum Q. In other systems, the deviations from the ideal boson features may influence the magneto-association 

C . of atoms into molecules via Feshbach resonances [3, llOl • 

^^ ■ Recently, several studies have focussed on the links between the composite particle nature of bosons and entan- 

t:]" . glement effects |ll| - [l3| using the principles of quantum information theory. It was shown that effects of the Pauli 
^D ' exclusion principle diminish when entanglement features dominates [ll|, [12, flil - llq . Measures such as purity P [12, Il7| 
^^ ] have been proposed to quantify the strength of entanglement between the constituent fermions. While the exclusion 
principle imposes the requirement that an antisymmetric state vector be assigned to an identical group of fermions, 
the rules becomes less stringent when the fermions become entangled and lose their fermion identity. The exact 
association between entanglement and the exclusion principle is not well understood as the latter implicates some 
degree of nonlocal interactions for fermionic particles of the same state to display anti-bunching behavior. Beyond 
C^ , a critical dimension, the tendency of similar fermions to "avoid" each other is lost, with the problem very much 
dependent on the system parameters such as dimensions and thermodynamic factors. Currently, very few details 
exists on the nature of the quantum correlation (whether classical or non-classical) that underpins the action of the 
exclusion principle in composite boson systems. 

In this work, we examine the conversion of composite bosons to orthogonal fermionic fragment states due to 
dominance of Pauli exclusion at increasing densities of the correlated boson system. The specific case of excitons is 
considered, and the ionic conversion of e — h — >■ X is examined as an equivalent of the well known excitonic Mott 
transition where the excitonic system ionizes due to space filling effects [13, [3 ^^ tiig^ densities, resulting in a free 
electron- hole plasma X. The situation in which the bound exciton state vanishes and merges into a continuum of 
scattered states can also be considered as the Mott transition. In this case, the scattered states exist as quasi-free 
fermions. In a recent study [201 , the Mott effect was noted to occur when the excitonic Bohr radius became nearly the 
same as the screening length, that is, the bound electron-hole pair states co-existed with the high-density electron- 
hole liquid phase. While the interplay of several many-body effects are responsible for the Mott transition, questions 
remain as to whether the conversion from the excitonic state to fermionic plasma-like phase proceeds in a continuous 
manner. 

Different timescales can be seen to operate in the composite boson-fermion system: the highly entangled elementary 



bosons which operate at fast times compared to the scattering times specific to fermionic species. To this end, the 
composite bosons can be analyzed via the open quantum system approach based on the system-plus-reservoir model. 
The total Hilbert space is divided into subsystems according to different timescales and/or states which exist in distinct 
phases. The Fock-Hilbert space of state vectors associated with a system of identical bosons (fermions) is spanned 
by only symmetric (antisymmetric) functions. The completeness theorem [2l| applies to many-particle state vectors 
in distinct Fock-Hilbert subspaces of antisymmetric and symmetric particle functions. However fluctuation in the 
number of interacting particles occur when energy dissipates from one subspaces to another. In this work, we consider 
that the Hilbert space of state-functions is decomposed into a composite boson subspace and an orthogonal fragment 
subspace of quasi-free fermions. The latter subspace enlarges as the composite boson constituents deviate from ideal 
boson commutation relations due to enhanced fermionic features. The quantum master equation approach of the 
Lindblad form [22|, [2j] may be used to investigate composite boson systems, however evaluation of the density matrix 
of such high dimension al sy stems will present insurmountable challenges. Stochastic formalisms such as the quantum 
trajectory method [24 . |25| involving non-Hermitian terms which cause quantum jumps, coupled with the Feshbach 
projection-operator partitioning technique |26| may provide a viable route to study the dynamics of composite bosons. 
In this study, we employ the Green's function formalism J27ll28l| to examine the dynamics of composite boson systems. 
A coupled system of boson condensate and fragment states with variable fermionic character has implications 
in the field of quantum information and processing. Composite bosons can be studied in the context of quantum 
tunneling of macroscopically coherent systems such as the double-well Bose-Einstein condensate. The double-well 
condensate is a well-known lattice system capable of variable controls 29|j_ 30| , exhibiting a range of quantum effects 



such as self-trapping, Josephson oscillations [3l| , and entanglement [32|, [3j] . In this work, we report on oscillations 
that occur in the coupled composite boson system, with testable predictions based on the quantum dot excitonic 
systems. These results are extended to analyze multichromophoric macromoleculcs (MCMMs) in photosynthetic light 
harvesting complexes. MCMMs are systems of great interest |34| - [38| due to the appearance of long-lived quantum 
coherences, even at physiological temperatures. 

This paper is organized as follows. In Section [11] we provide a brief review of cobosons or composite boson states 
and their Schmidt decomposition properties, along with description of the Schmidt number and purity measures. The 
occurrence of the fermionic fragment state which lies orthogonal to the composite boson is highlighted. In Section Hill 
a description of other alternative measures of deviations from ideal boson characteristics is provided, and numerical 
estimates of the bosonic deviation measure in quantum dots of varying size and fermion pair number is provided. 
In Section IIVI we introduce the open quantum system model and master equation for the cxciton-fermion system. 
The importance of Non-Markovian dynamics due to the fermionic background is briefly described in this Section. 
In Section |Vl the tensor structure of the many particle Fock-Hilbert space is examined for the N excitonic bosons, 
and the tunneling dynamics of composite bosons states is numerically examined using the decay branching ratio in 
Section rvTl A discussion of the appearance of exceptional points under experimentally testable conditions is included 
as well. In Section [VII[ we examine the application of results obtained in Section IVll to pigment protein complexes in 
light- harvesting systems, and present our conclusions in Section IVIIII 

II. SCHMIDT DECOMPOSITION OF COBOSONS OR COMPOSITE BOSON STATES 

Following earlier formalisms of cobosons comprising two fermions [ll|, [l2|, [S^] , the coboson creation operator of 
distinguishable fermions in the Schmidt decomposition involving a single index is written as 

5'=EVA7a]6], (1) 

where Xj are the Schmidt coefficients, a, and b] are different fermion creation operators in the Schmidt mode j, and 
J denotes the number of Schmidt coefficients |40| . The distribution of Xj is linked to a measure of entanglement, 
which is provided via the Schmidt number |ll| 

CJO 

IC^l/Y^X^ (2) 

The Schmidt number K. is the quantum counterpart to the classical Pearson correlation coefficient, and is an important 
entanglement measure [4l|,|42] where large K. indicates high correlations and entanglement. K. is also linked to another 
quantity known as the purity P [l2| via ^=-f, and varies between zero and one. In the case of two particles, P = Tr p^, 
where p is the density matrix of the examined particle. 



For two identical particles, (fermions or bosons), the symmetrization postulate constraints a boson (fermion) state 
associated with the system to be totally symmetric (antisymmetric) under permutation of the particles. As a con- 
sequence, the Schmidt decomposition of the state involves more than one term, and a bipartite state of two indis- 
tinguishable particles is generally considered entangled. This highlights the importance of the Schmidt number in 
determining entanglement in quantum states of identical particles [43|. B and B^ obey the non-bosonic commu- 
tation relation, [5,5^] = 1 + sA, where s = +1(— 1) if the two interacting particles are bosons (fermions), and 

A = ^ ^j [ olflj + bjbj ] . The state of N composite bosons appear as 
j=o ^ ' 

|7V) = ^^i=-|0), (3) 

where deviations from the ideal boson characteristics arc contained in the normalization term xn obtained using 
(A^|A^)=1. The effectiveness of the operator i? as a bosonic annihilation operator can be seen via the action of 
operator B on state |A^) [il| 



B|Ar) =aArVAr|A^-l) + |J-^), (4) 

where ajq is a ideality parameter and \Tjss) is the fragment state resulting from the non-ideal nature of |A^). aj^ = 

. ^" , which in the case of an ideal composite boson yields the normalization ratio Xn±i/xn -^ 1- This ratio is 

seen as a measure of the degree of ideal bosonization for a state of N cobosons, and appears in the pair number mean 
value corresponding to the state |A^) as 

W . <^.l+(„-l,«±l (5) 

{N\N) XN 



and also in the commutator mean [ll|. 



{N\[B,B^]\N)^2^^^^^-l (6) 

XN 



A neat inequality involving the upper and a lower bound to the normalization ratio was determined as 1 — P • A'^ < 

^^^^ <l-P [H. 

The fragment states, \J-n) remain orthogonal to |A^ — 1), yielding the correction factor [11U39J 

(j-^l^^) = i-^^E±l^N (^^ - ^^"j . (7) 

XN \XN-1 XN J 

The ratio Xn+i/xn is strictly no n- increasing as A^ increases |12| . and for small bosonic deviations such that 
Xn+i/xn ~ 1 — ^, the last term in Eq. ([7]) can be dropped, {J-n\^n)~ S and the commutator mean, {N\[B,B'^]\N) = 
1 — 25. In the limits, ajv — >■ 1, 5 — >■ 0, {J-n\^n) —>■ 0. The formation of the fermionic fragment can be compared to the 
formation of an electron-hole plasma when the density of a collection of correlated electron-hole is increased, giving rise 
to an enhancement in fermionic features. With increasing closeness of interacting paired fermions, the electron-holes 
pairs become unbound as is the case when lattice temperature is also increased. The state of A^ composite bosons (see 
Eqs. ([3])) with a well-defined atom number evolves into a mixture of lower number states and fragment state \J^n), 
characterized by the fidelity decay, 7. The role of 7 is significant, due to its influence in a Zeno-like mechanism where 
repeated measurements halts further decay of the composite boson states. 

III. ALTERNATIVE MEASURES OF DEVIATIONS FROM IDEAL BOSON CHARACTERISTICS 

Here we consider other measures that can be used in place of the purity P. An experimentally accessible measure 
that can be used to capture deviations from ideal boson characteristics is based on the (normalised) second order 
correlator (72 H: EH , which characterizes the probability of detecting of particles at times t and t + t 

^ {BHt)BHt + T)B{t + T)B{t)) 
{N{t)){N{t + r)) 



g2{T) is based on correlations of the boson operators where N is the number operator, N\n) 
N ferniion pairs. At zero delay, 



52(0) 



(JVo(t)(iVo(t) - 1)) 

(iV)2 



n\N) associated with 



(9) 



where {No{t)) ^{B'^ {t)B{t)) . The seeond-order correlator provides information on the underlying statistical features, 
such as the Poissonian case (52(0) = 1) in coherent systems involving a large number of Fock states. The anti- 
bunching , sub-Poissonian case (172(0) < 1) is applicable in the fermionic limit at high densities. The classically 
accessible thermal states which display a bunched, super-Poissonian distribution gives rise to (72(0) > 1, we do not 
consider such states in the work here. A simple form of the zero-delay t ~ Q correlations for the single-mode state, 

was obtained as 32(0) — ""^^2"" Hi where ajv = ViV^l — 2(A^ — 1)(^)2. The latter expression is applicable to 
excitons of bohr radius a-Q in quantum dots of size L at small values of ^ and n <C — . In the pure bosonic case, 
<72(0) = {N — 1)/N -^ 1 as N -^ 00. The fermionic structure of excitons becomes noticeable with decrease in 32(0), 
here we estimate the bosonic deviations using 5 = 1 — 32(0). Results displayed on Fig. [T]show the increase in the 
bosonic deviation measure 5 with increase in number density (or decrease in quantum dot size) as quantified by the 
^. These results translates to the growth of overlap in fragment states, {Fn\^n)^S with increase in ^. 



ratio, -^ 




A^ 



FIG. 1: Bosonic deviation measure 5 as function of fermion pair number A'^ in quantum dots at varying ratios, — 
solid line), 0.03 (Blue), 0.05 (Red), 0.07 (Green) 
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An alternative measure of the bosonic deviations is obtained using the degree of binding of the paired fermions 



ad 



Ec 

1 

Eb 



(10) 



where Eb is the maximized binding energy of the composite boson, and Ec is the smaller binding energy of the quasi- 
bound fermion pairs, which reaches zero for free fermions. This definition is intimately linked to the distinguishability 
of paired fermions, a system of many strongly bound fermion pair is less distinguishable and more entangled than 
one of free fermion species. In the case of excitons in semiconductors at varying temperatures and densities, a mixed 
phase of bound excitons and free electronhole plasma is formed. Here the degree of ionization of the electronhole 
plasma, a; = 1 — ■^, is a suitable candidate for gauging the bosonic deviations of the excitonic system. The term 
nb denotes the excitonic density and Uf denotes the density of the free electrons or holes at a given temperature and 
density. The evaluation of a^ requires carrier density parameters such as the chemical potentials of distinct species, 
temperature, band multiplicity, spin degeneracy and an integrand with the retarded Greens function of carriers [20[. 
Yet another quantity that provides a measure of deviations from ideal boson characteristics is the iV-particle non- 
escape probability in which N paired fermions are found within the same region, E 



PN{t) 



N 

/ J|dr„|^(ri,...,rw;t)|^ 

'J z2 „ 1 



(11) 



The evolution dynamics of PN(t) can be examined further by considering the tensor structure of -PAr(i) in Fock-Hilbert 
space, however this procedure involves the incorporation of all degrees of freedom of the A'^-particle system. In Section 
IVl we employ a similar approach, but one which requires fewer parameters, applied to the case of the excitonic bosons. 



IV. OPEN QUANTUM SYSTEM MODEL AND MASTER EQUATION FOR EXCITON-FERMION 

SYSTEM 

We first examine tlic case of a paired electron-hole or exciton interacting with a background of dissociated electrons 
and holes. The total Hamiltonian of the exciton and fermion reservoir appear as 

-Ht = -He + nf + m + nd, (12) 

-He = J2^E,-ijb)b],Bk (13) 

K 



Uf = ^(ee - He)a\^ ttk, + ^(eh - Mft)/»L ^/^i^ ^ '^'s XI alhlk'^-k'ttk' (14) 

kc fch A;, A;' 

U, = ^a*(fce,fch)aL^-fc,^..-., +'^(^c,fch)i?I_,_^afe„/»-fe, (15) 

kc ,ki^ 

-Hd - z5^7^(fce,fch)BV.,^.e-.. (16) 

where the subscripts, e, / in the Hamiltonian terms refer, respectively to the exciton and fermionic fragments, T-Li 
(Eq. P^ ) denotes the interaction between the two subsystems, aj. (h\) denote the electron (hole) creation operator 
with wavevector A:c(fch) and kinetic energy ee(eh)- The boson creation operator, Bj^ is labeled by the wavcvector K 
and is obtained using the Fourier transform of the site-dependent operator, Bj as B*j^ = N^^^'^'^^e^^^-^ Bj . The exciton 
creation operator Bj^ which is localized in i^-space, is delocalized in real space. Ec in Eq. P^ may be seen as the 
minimum energy required to form the composite boson system, /ig and fih denote the respective chemical potentials 
for the electrons and holes, while fin denotes the chemical potential of the bosonic system. 7^ denotes the strength 
of electron-hole interactions which occurs in the fragment subspace. We exclude, for simplification purposes, the 
background plasma of electron-hole carriers that are originally formed when excitons are created, a in the interaction 
operator Hi (Eq. (fT5|) ) is the momentum dependent coupling strength between two dissociated fermions and the 
correlated fermion pair that constitutes the exciton. The non-Hcrmitian dissipative Hamiltonian Hd (Eq. (|16p ) is 
dependent on the spontaneous decay jn, which increases with greater deviations from ideal bosonic features. This 
term may be attributed to the growth of the free fermion plasma state, or fragment state (see Eq. (O), resulting in 
the breakdown of boson states. 

The density operator p of the quantum system associated with the total Hamiltonian, Ht (Eq. ([T2|) ) is obtained 
from the generalized Liouvillc-von Neumann equation -^ = —iCp. The generator C ^ Cb + Ci + Cf maps the initial to 
final density operators via a Liouvillc superoperator $(t, 0): p{0) 1— > p{t) = $(t,0]p(0). The Liouville-von Neumann 
equation can be recast as a master equation of the following Lindblad form [221 . |23| : 

^p{t) = -tint + n.,p{t)] + J2 Ik (vkpml - kv^v,, p{t)}\ , (17) 



fc=l ^ 

where Vk denote Lindblad operators, in which both an operator and its hermitian conjugate are labeled by k, and 
the decay terms, {7*:} constitute the spectrum of the positive definite d-dimensional Hermitian Gorini-Kossakowski- 
Sudarshan matrix A |23| . The first term in Eq. (J17p represents reversibility in system dynamics, while the sym- 
metrized Lindblad operators, Vk involve transitions between the many-particle levels of both the composite boson 
and states present in the fermionic background. While the notations associated with these possible energy transfer 
processes are not shown in Eq. (1171) . we note that these processes contribute to a range of dynamical time scales 
due to the possible number of energetic degrees of freedom that can arise. The multitude of these transitions adds 
to the complexity of solving many coupled differential equations of Eq. (|17p . The problem is however tractable in 
systems with weak system-reservoir coupling when the Markov approximation applies. Next, we briefiy describe the 
non-Markovian dynamics that may occur in the composite boson system. 

A. Non-Markovian dynamics due to the fermionic background 

At the initial period of quantum evolution, the dynamics of the composite boson-fermion system is dominantly 
non-Markovian, determined by the fermion bath memory time, and the Lindblad form in Eq. (J17p does not provide 
a reliable description of the dynamics. Several processes may give rise to the non-Markovian dynamics, including 
the dynamics associated with the vibrational environment of the composite boson-fermion system. The time-scales 



of processes which occur in ideal and composite bosons differ due to decreased scattering between two elementary 
bosons as compared their composite counterparts |46| . There also exist differences between ideal bosons and composite 
bosons in terms of the decreased scattering via lattice vibrations in non-ideal boson systems due to phase filling of 
the fermionic phase-space. 

Information about the short time dynamics of the composite boson-fermion system will help in the understanding 
of the subtle links between non-Markovian dynamics and entanglement measures such as the Schmidt number /C and 
purity P of the parent composite boson state. This may contribute to accurate calculations of binding energies of 
complexes consisting of several fermions such as excitonic complexes (47| . and improvements in density functional 
approaches to determining electronic properties of material systems. There are inherent difficulties in a detailed 
analysis of non-Markovian features, as this will require the use of a non-Lindblad set of relations which incorporate 
the finite time scale of the vibrational modes of the uncorrelated fermions. To keep the problem tractable, we consider 
in the next Section, the conversion of the composite boson state \N) to the fermionic fragment state \J^n) using a 
simple open quantum system consisting of a two-level system interacting linearly with a dissociated electron-hole 
(fermion) reservoir. We note that as the density (or lattice temperature) increases, the forward conversion of boson 
into free fermions is favored due to screening of the Coulomb interaction that tends to bind and form composite 
bosons. Hence, at higher densities, the interaction between the boson-fermion background system is more likely to be 
Markovian. 

V. TENSOR STRUCTURE OF THE MANY PARTICLE FOCK-HILBERT SPACE: N EXCITONIC 

BOSONS 

The tensor structure of the many particle Fock-Hilbert space for the system of N correlated electron-hole or 
composite exciton state appear as 

ST=Sb®Sf (18) 

where the total many particle Fock-Hilbert space St is expressed as the tensor product of two orthogonal subspaces. 
Here S^, {Sf) is the subspace associated with the composite boson state |A^) (fragment state, {Tn) ). Both subspaces 
are considered to include the complete set of bound and unbound states, and combination thereof to incorporate 
interactions present in any A^-body system. 

To keep the problem tractable, we assume that all N bosons possess the same wavevector, K — Kq, and consider 
that the subspace Sb is spanned by (|l)b, |0)fe), where |l)fe denotes the presence of A^ composite bosons with wavevector, 
Kq, and |0)b indicates a lower {N' < N) number of bosons within Sb- As noted in Section ITVl the Fock space of 
composite boson and fermions allows for many complicated interactions, coupled with the inseparability of the degrees 
of freedom of the boson system and its fermionic background. The various correlations that occur between |A"') and 
I A^") (A^', A^" < A^) is implicit in the defined jO);,, and hence we consider a collective state that involves a superposition 
of similar states, 

|0), = J2 MN')\N') (19) 

N'<N 

A{N') is a weight factor that is dependent on the evolution dynamics of the coupled boson-fermion system, and for 
simplification we ignore details of the possible interactions between the different composite boson states. The raising 
and lowering operators appear respectively, cr^ =|l)hf,(0| and cr^L =|0)f,f,(l|. Due to the choice of definitions, a crude 
estimate of the energy difference between |l)f, and |0)h states is given by the binding energy of a single composite 
boson. This quantity can be obtained using experimental techniques in the case of excitonic systems. 

Likewise, we also consider that Sf is spanned by (|1)/, |0)/) where |1)/ (|0)/) denotes the presence (absence) of 
the fragment state, \J^n)- Instead of identifying electron and hole states, we denote both types of fermions using 
operators cl,,Ck, and the raising fcrnfionic operator, cr{^_ appear as linear combination of cj. operators 

'^i = |l)//(0| = ^E-?'4 (20) 

which can be easily shown to obey the anti-commutator relation {a_,cr!^_} =1, where the lowering operator, cr_ = 
|0)//(1|. ^=X]fe 7fc' where jk is the weight amplitude for the respective fermion in simplified system of noninteracting 
fermions. 

Using Feshbach projection-operator partitioning method [2g, the total Hilbert space of "Hr (Eq. ((T^ ) is divided 
into two orthogonal subspaces, Sb and Sf (see Eq. p^ ) generated respectively by a projection operator, V—a^d'L 



and its complementary, Q= 



J^f 



with Hence PQ=QV=Q, and the reduced density state operator associated with 



the central composite boson system of interest is obtained via 

Pb = V ptV, 
= TvfipTit)} 



(21) 
(22) 



where pT is the density operator of the total system described by Ht- The reduced density state, pb can also be 
obtained by taking the partial trace over the fermionic environment as shown in Eq. (|22p. The density operator 
associated with the fermionic fragment can be obtained using Pf{t) ~ Q prQ- If we consider that at time to, the 
subsystems of composite bosons and fermionic fragment are in separable states, the evolution oi pi,{t) in Laplace space 
becomes 



Pb{z)-[z-{Cb + N{z))] 



(23) 



where the energy term N(z) is associated with non-Markovian interactions. Feshbach projection-operator partitioning 
has been employed in an earlier work [48|, to examine the dynamics of open quantum systems via the stochastic 
quantum trajectory approach. While the stochastic route presents a physical interpretation of the quantum trajectories 
in the case of Markovian dynamics, it offers no viable explanation for the occurrence of non-Markovian dynamics, 
due to the finite correlation time of the non-Markovian reservoir [49| . An alternative approach involving the post- 
Markovian master equation j49l |50| is known to be applicable in the regime between Markovian and non-Markovian 
quantum dynamics in open quantum systems. Here we utilize the Greens function approach [2 71 128| to provide an 
effective description of the quantum evolution of the composite boson system. 

VI. TUNNELING DYNAMICS OF COMPOSITE BOSONS STATES: APPEARANCE OF 

EXCEPTIONAL POINTS 



We consider the tunneling dynamics between a pair of A^-composite boson states |l)6i and \l)b2 with a total 
Hamiltonian of the form (^=1) 






^b2h2 
■ ^b2 cr+ cr_ 



T/* „blh2 



Va'Ja 



b2hl 



*7di<'^" 



^Jd2<7^'<J'' 



(24) 



where uJbi (i=l,2) are the two composite boson transition energies (c.f. Eq. (fT3|)). V denotes the tunneling energy 
between the two iV-composite boson states, and is taken to be real and positive, without loss in generality. Due to the 
choice of definitions for |l)b and |0)b, the tunneling energy V could involve the transfer of excitation associated with 
just one fermion pair. As will be described in Section IVHl V could represent a series of repeated processes leading 
to transfer of states from one site to another in photosynthetic protein complexes. The dissipative terms in Eq. (|24|) 
represent leakages of boson states into the fermionic subspaces, Sb —> Sf. The state |l)(,i(|l)(,2)decays at the rate 
7di (7^2) to the lower state |0)6i (|0)62)- The phenomenological rate jdi=^i Si, is employed where Si={TN\J^N) (see 
Eq. ([7])) is a measure of the deviation from bosonic features with increase in density of the composite boson state in 
a specific subspace Sbi- A^ is taken as a constant with units of energy. 
We consider the retarded Green's function of the form 



where 9(i) denotes a step function. The Fourier transform, Gi^2{E)=J 
appear with the Lippmann-Schwinger matrix terms 27 1 



Gi.2(t)--*e(t)({a^i(t),af(0)}) (25) 

dtGi,2(i)e*^* for the system in Eq. ((24|) 



GVliE) 



E - uJbi + iri+ -V 

— V E — ujb2 + ii]- 



l( Idi 
2 V 7d2 



(26) 



where 77 is a very small number and the dissipative process associated with 7^1 (7^2) in the subspace Sbi {Sb2) is 
considered as an irreversible loss of the composite exciton of size N . In an initial state, i = and the fermionic 
fragment will be absent, and only the N-stsXe composite boson state |l)hi is excited. We denote the probability of 
excitation to remain at its initial site, Pi.i(0)=l. The probability Pi_2 of excitation transfer from one boson subspace 
to another, Sbi -^ Sb2 is determined by inverting Eq. 



2^2 



PiMt) = 77^e-'^^*(coshfi,f - cosf},.t), 



(27) 



where 7d=|(7di+7d2), ^ = i^r + i^i = \/^V^TJto^^^vydf , wo=Wb2-w&i,7d= ^(7^2-7^1)- In the absence of the Pauh 
exclusion related dissipation at time t > t' , the boson- fermion system undergoes Rabi-type oscillation determined by 
Wo and V. In the continued presence of dissipativc terms, the total probabilities, P1.1+P1.2 < 1 is not conserved and 
there is loss of normalization which is dependent on 7^1,7^2- 

Fig. [2ji,b shows the tunneling dynamics for the specific case when the energy difference a;o=0. Depending on 
the tunneling energy V and decay rates F', there is existence of a coherent regime {2V > jd) or incoherent regime 
{2V < 7rf). At large V and small dissipation levels, there is increased exchanges between the two coupled bosonic 
states (Fig. [5^) . The gradual decrease in the A^-state boson population appears due to increased deviation from ideal 
boson characteristics associated with increase in number density (or decrease in quantum dot size) as shown in Fig.[2]3. 

The appearance of topological defects known as exceptional points [5l| occurs when 17 = 0, V = ^ Unlike 
degenerate points, only a single eigenfunction exists at the exceptional point due to the merging of two eigenvalues. 
The critical boson densities and lattice temperatures at which exceptional points occur is evaluated using many body 
theory which takes into account dynamical screening processes. This screening may arise from Coulomb interactions of 
the one-particle and two-particle properties between the same and different fermion species constituting the composite 
boson system. These special points may be associated with a range of system parameter attributes, and hence lie 
within an allowed spectrum that may be amenable to experimental detection. 
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FIG. 2: (a)Population difference, AP—Pi,i-Pi.2 as a function of time t, and coupling energy, V at dissipation rates 7^1 = 
7ii2=0.1, for the degenerate case (a;62='^6i). The units are chosen such that h=l, Ai=A2 = 1 (i.e 5i=S2=0.1). Time t is 
obtained as multiple of to, the inverse of Qo (at 7^1 = 7d2=0.1) (b) Population difference, AP—Pj^i 
t, and deviation factor, S = S2 at tunneling energy V=l, 7di=0 and A2 = l. 



Pi, 2 as a function of time 



The decay branching ratio is quantified by the fraction Fi (or F2) of a A^-state composite boson that decay via jdi 
(or 7^2)- F2 is evaluated using Parseval's theorem [52] 
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The fraction Fi = l — ^2- Results displayed on Fig. [5]i show the increase (decrease) in the branching fraction F2 due to 
increase (increase) in the bosonic deviation measure 62 (Si) for fixed values of the energy difference ujq and tunneling 
energy V. Fig. [3}d shows the notable decrease of the branching fraction F2 with increase in the energy difference luq. 
Conversely, the branching fraction Fi increases with increase in the energy difference ujq, as expected. 

VII. APPLICATION TO PIGMENT PROTEIN COMPLEXES IN LIGHT-HARVESTING SYSTEMS 



It is useful to analyze the results obtained in Section I VII in the context of large photosynthetic membranes which 
constitute many biological pigment-protein complexes (i.e., chromophorcs) such as FMO (Fenna-Matthews-Olson) 
complexes in the green sulphur bacteria [3J, |35| . The FMO complex trimcr is made up of three symmetry equivalent 
monomer subunits, with each unit constituting eight bacteriochlorophyll (BChl)a molecules supported by a cage of 
protein molecules. The FMO complex acts as an efficient channel of excitation transfer in which photons captured in 
the chlorosome which is the main light harvesting antenna complex, are directed via a series of excitonic exchanges to a 
reaction center (RC) where energy conversion into a chemical form occurs. Long-lived coherences between electronic 
states lasting several picoseconds, much shorter than the 1 ns dissipativc lifetimes of excitons |35| . have been a 
topic of intense investigation in recent years [36-38;. i53l | . Currently it is still not clear as to how biological systems 
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FIG. 3: (a) Branching fraction F2 of a TV-state composite boson that decay via 7^2 as a function of bosonic deviations, 5i 
and 62- Energy difference uo=uji,2-ijJbi=0.5, and tunneling energy V=l. (b) Fraction F2 as a function of energy difference uq 
and bosonic deviation, 62- The tunneling energy V=5 and (5i=0.1. The units are chosen such that h=l, Ai=A2 = 1 for both 
figures. 



comprising hundreds of photosynthetic complexes and many more correlated excitonic states act in unison to maintain 
the quantum coherences in the noisy environment, and attain the much envied high efficiencies at psychological 
temperatures. 

In the FMO complex, the chromophore sites numbered 3 and 4 are located near the reaction center, and thus are 
closely linked to the sink region where energy is released, while chromophore sites 1 and 2 are strongly coupled, and 
dissipate energy via site 3 |3J]. The sites 1, 6, and 8 are located at the baseplate which connects to the chlorosomes 
that receive electronic excitation. Recently, it was shown that the eighth chromophore (at site 8) is located nearer the 
chromophore sites of neighbouring monomers compared to sites in its own monomer [53| • This indicates a stronger 
inter-monomer (as compared to intra-monomer) interaction as far as the eighth chromophore is concerned and the 
excitation at site 8 is likely to propagate to a inter-site monomer (see Fig. 2^ ). It was pointed out [5J| that the eigth 
chromophore acts to facilitate excitation transfer between monomers of the FMO trimer even though it is best placed 
to receive excitation at the earliest time. To this end, the eighth chromophore appears to play a critical role in the 
topological connectivity of large molecular structures of multichromophoric macromolecule (MCMM) systems. The 
electro-optical properties of these MCMMs vary from those of single chromophores, depending on the delocalization 
of excitons within each MCMM. 

Here we apply Eq. ((26|) to a model in which excitons in MCMMs undergo tunneling dynamics over distances that 
are large compared to the average distance between chromophores within the FMO monomer/trimer configurations 
[56j . A simple setup is shown in Fig.HjD, where 1 < A^ < Nm—'iA, and which describes the system of composite bosons 
in one FMO trimer that is considered as the multichromophoric macromolecule. Individual MCMM can be coupled 
to each other as is shown for two FMO trimers in Fig. [JJd. Dissipation may arise from Pauli exchanges between 
the FMO pigments and the protein bath, and recombination and trapping effects specific to each macromolecular 
system. The configuration Fig. Hb can be further extended to one in which each MCMM includes several trimers 
forming aggregates. The net dissipation within each MCMM may vary from other similarly configured MCMM, 
depending on the connectivity and proximity to the region of light illumination. In a recent work ^55|], multiple 
detrapping/retrapping processes as opposed to the slow (200 ps) direct transfer between RCs in the purple bacterium 
Rhodohacter sphaeroides [55| were noted to contribute to the delocalization of excitation among several reaction 
centres (RCs). In this regard, the tunneling energy V in Eq. (j26p may be based on a cumulative process of repeated 
trapping/detrapping events instead of a single direct transfer mechanism. 

Oscillations between MCMMs are expected as shown in Fig. [2^,b, with excitation exchanges that gradually fades 
with time depending on the initial conditions, dissipation parameters, 7^1,7^2 (specific to the two MCMM sites) 
and average energy difference (« ujq) between the MCMM sites. As noted earlier, coherence times are much shorter 
than the dissipative lifetimes of excitons [35|. Of particular interest is use of the branching ratio in Eq. (j29p . which 
identifies effective routes of energy propagation in large topologically connected network structures. A single antenna 
complex may serve several FMO complexes, and a reaction center may be linked to several FMO trimers. Selective 
MCMM sites which function as collection centers may experience greater dissipation than other sites and possess 
higher branching ratios (Fig. [3^,b). Such sites contribute to a structural arrangement that enable photosynthetic 
organisms to better utilize cellular resources. 

Continued coherence in photosynthetic complexes is ensured by a small bosonic deviation measure S and large 
number A^ of fermion pairs involved during energy exchanges (Fig. [2]). This suggest that a molecular environment 
that is highly correlated, with a large Schmidt number /C (Eq. 1^ is likely to preserve electronic coherences needed 
for propagation of excitation such that energy is harvested efficiently. The spectral and molecular dynamics of 
MCMM at various lattice temperatures and excitation densities (light illumination) appear to influence the kinetics 
of exciton and electron-hole pair recombination and relaxation processes. These can be seen as critical factors that 
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influence the Schmidt number K. specific to a MCMM. Thus far, we have considered a few factors which underpins 
the high efficiencies noted in light-harvesting systems, on qualitative terms. A quantitative approach would involve 
modeling the realistic condition of an entire photosynthetic membrane constituting many FMO complexes, and taking 
into account the topological connectivity of thousands of bacteriochorophylls. The approach taken in this work is 
expected to help understand the importance of the coexistence af boson and fermionic phases, Pauli scattering effects 
and selective dissipation in photosynthetic systems. 
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FIG. 4: (a) Simplified trimer configuration in which the eighth chromophore is positioned close to the first chromophore of a 
neighbouring monomer, resulting in a stronger inter-monomer (as compared to intra-monomer) interaction for an excitation at 
site 8 [5J] . (b) A system of two coupled FMO trimer multichromophoric macromolecules that is capable of displaying long-lived 
coherences for small bosonic deviations of composite excitons at each trimer site. 



VIII. CONCLUSION 

In conclusion, the quantum dynamics of conversion of composite bosons into fermionic fragment species is demon- 
strated using an open quantum system approach based on a system-plus-reservoir model. The total Hilbert space 
which constitutes a composite boson subspace and an orthogonal fragment subspace of fermions pairs is used to ex- 
amine the effect of system parameters during the tunneling dynamics of coupled composite bosons states. The results 
highlight the interplay of boson and fermionic phases that is dependent on density (and indirectly the network connec- 
tivity) and lattice temperature, and the appearance of exceptional points based on experimentally testable conditions 
(densities, lattice temperatures). The effect of Pauli exclusion and other dissipative factors on the multichromophoric 
macromolecules (MCMMs) in photosynthetic light harvesting is examined in the light of quantitative results obtained 
for coupled composite bosons systems. It is noted that long-lived quantum coherence in photosynthetic complexes is 
assisted by small bosonic deviation measures (large Schmidt number) and large number of excitons involved during 
energy exchanges, to give rise to a highly correlated molecular environment. Moreover, specific MCMM sites which 
function as collection centers may possess higher branching ratios, and contribute to a structural arrangement that 
enable photosynthetic organisms to better harness solar energy efficiently. 
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